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Abstract 

We investigate three-colour QCD thermodynamics at finite quark chemical 
potential. Lattice QCD results are compared with a generalized Nambu Jona- 
Lasinio model in which quarks couple simultaneously to the chiral condensate 
and to a background temporal gauge field representing Polyakov loop dynam- 
ics. This so-called PNJL model thus includes features of both deconfinement 
and chiral symmetry restoration. The parameters of the Polyakov loop effec- 
tive potential are fixed in the pure gauge sector. The chiral condensate and 
the Polyakov loop as functions of temperature and quark chemical potential 
are calculated by minimizing the thermodynamic potential of the system. The 
resulting equation of state, (scaled) pressure difference and quark number den- 
sity at finite quark chemical potential are then confronted with corresponding 
Lattice QCD data. 
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1 Introduction 



Recent years have seen an expansion of activities devoted to the study of the QCD 
phase diagram. Heavy-ion experiments are looking for signals of the Quark-Gluon 
Plasma. Large-scale lattice simulations at finite temperature have become a princi- 
pal tool for investigating the pattern of phases in QCD. Accurate computations of 
lattice QCD thermodynamics in the pure gauge sector have been performed. First 
results at finite quark chemical potential are available. The equation of state of 
strongly interacting matter is now at hand as a function of temperature T and in a 
limited range of quark chemical potential //. Improved multi-parameter re-weighting 
techniques [1,2], Taylor series expansion methods [3-5] and analytic continuation 
from imaginary chemical potential [6-9] provide lattice data for the pressure, entropy 
density, quark density and selected susceptibilities. 

A straightforward interpretation of these data in terms of QCD perturbation 
theory does not work because of poor convergence at any temperature of practical 
interest [10, 11]. In order to overcome this problem, resummation schemes have been 
proposed, based for example on the Hard Thermal Loop (HTL) approach [12-20] or 
on dimensionally reduced screened perturbation theory (DRSPT) [21-23]. However, 
these approaches still give reliable results only for temperatures T > 2.5 T c , far 
above the critical temperature T c ~ 0.2 GeV. At these high temperatures, the HTL 
approach motivates and justifies a picture of weakly interacting quasiparticles, as 
determined by the HTL propagators. 

In order to extend such descriptions to lower temperatures closer to T c , various 
models have been proposed. Early attempts were based on the MIT bag model [24]. 
More sophisticated approaches became necessary when more precise lattice data ap- 
peared. Various aspects of QCD thermodynamics have been investigated in terms 
of quasiparticle models based on perturbative calculations carried out in the HTL 
scheme [25-30], in terms of a condensate of Z 3 Wilson lines [31], by refined quasi- 
particle models based on the HTL-resummed entropy and extensions thereof [32], 
by an improved version with a temperature-dependent number of active degrees of 
freedom [33,34], by an evaporation model of the gluon condensate [35], by quasi- 
particle models formulated in dynamical terms [36], and by hadron resonance gas 
models below the critical temperature [37] (for a recent review see [38]). 

In this paper, we study the thermodynamics of two-flavour QCD at finite quark 
chemical potential. Our investigation is based on a synthesis of a Nambu Jona- 
Lasinio (NJL) model [39-44] and the non-linear dynamics involving the Polyakov 
loop [45-49]. In this Polyakov-loop-extended (PNJL) model, quarks develop quasi- 
particle masses by propagating in the chiral condensate, while they couple at the 
same time to a homogeneous background (temporal) gauge field representing Polyakov 
loop dynamics. 

The "classic" NJL model incorporates the chiral symmetry of two-flavour QCD 
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and its spontaneous breakdown at T < T c . Gluonic degrees of freedom are "in- 
tegrated out" and replaced by a local four-point interaction of quark colour cur- 
rents. Subsequent Fierz transformations project this interaction into various quark- 
antiquark and diquark channels. The colour singlet qq modes of lowest mass are 
identified with the lightest mesons. Pions properly emerge as Goldstone bosons at 
T < T c . However, the local SU(N C ) gauge invariance of QCD is now replaced by 
a global SU (N c ) symmetry in the NJL model, so that the confinement property is 
lost. Consequently, standard NJL-type models are bound to fail in attempts to de- 
scribe N c = 3 thermodynamics around T c (and beyond) for non-zero quark chemical 
potential /i. 

On the other hand the NJL model, with just the simplest possible one-parameter 
colour current-current interaction between quarks, is remarkably successful in repro- 
ducing the thermodynamics of N c = 2 Lattice QCD at finite fi [50]. Encouraged 
by this result, the NJL quasiparticle concept does suggest itself as a useful starting 
point. However, whereas aspects of deconfinement are less significant in the N c = 2 
case, they figure prominently for N c = 3. This motivates our extension towards 
the PNJL Lagrangian as a minimal approach incorporating both chiral symmetry 
restoration and deconfinement. 

The deconfinement phase transition is well defined in the heavy-quark limit, 
where the Polyakov loop serves as an order parameter. This phase transition is 
characterized by the spontaneous breaking of the Z(2>) center symmetry of QCD 
[51-54]. In the presence of dynamical quarks the center symmetry is explicitly 
broken. No order parameter is established for the deconfinement transition in this 
case [55], but the Polyakov loop still serves as an indicator of a rapid crossover 
towards deconfinement. The chiral phase transition, on the other hand, has a well- 
defined order parameter in the chiral limit of massless quarks: the chiral (or quark) 
condensate (qq). This condensate, and its dynamical generation, is the basic element 
of the original NJL model. 

The primary aim of this paper is to test the effectiveness of the PNJL approach 
when confronted with Lattice QCD thermodynamics. The PNJL Lagrangian is 
derived in Section 2. Parameters are fixed in Section 3 by reproducing known 
properties of the pion and of the QCD vacuum in the hadronic phase, while the 
Polyakov loop effective potential is adjusted to pure gauge lattice results. Sections 4 
and 5 deal with the thermodynamics derived from the PNJL Lagrangian in the mean- 
field approximation. This discussion includes a detailed comparison with Lattice 
QCD results at zero and at finite chemical potential. 
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2 The PNJL model 



Following [47] we introduce a generalized Nf = 2 Nambu Jona-Lasinio Lagrangian 
with quarks coupled to a (spatially constant) temporal background gauge field rep- 
resenting Polyakov loop dynamics (the PNJL model): 

Cpnjl = $ (i^D* - m ) V + % [(W 2 + (^75^) 2 ] - W ($L4], Q[A], T) , (1) 

where = (ifj u , ijjA is the quark field, 

D^ = d"-iA^ and A" = S^A . (2) 
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The gauge coupling <? is conveniently absorbed in the definition of A^(x) = gA%(x 
where A% is the SU(3) gauge field and A a are the Gell-Mann matrices. The two- 
flavour current quark mass matrix is m = diag(m ui m d ) and we shall work in the 
isospin symmetric limit with m u = m& = itlq. A local, chirally symmetric scalar- 
pseudoscalar four-point interaction of the quark fields is introduced with an effective 
coupling strength G. 

The quantity U(<&, T) is the effective potential expressed in terms of the traced 
Polyakov loop 1 and its (charge) conjugate, 

$ = (Tr c L)/iV c , $ = (Tr c Lt)/iV c . (3) 

The Polyakov loop L is a matrix in colour space explicitly given by 

f/3 



L(x)=Vexp i [ drA 4 (x,r) , (4) 
. Jo 



with j3 = 1/T the inverse temperature and A4 = iA°. In a convenient gauge 
(the so-called Polyakov gauge), the Polyakov loop matrix can be given a diagonal 
representation [47]. 

The coupling between Polyakov loop and quarks is uniquely determined by the 
covariant derivative in the PNJL Lagrangian (1). Note that in the chiral limit 
(rho — > 0), this Lagrangian is invariant under the chiral flavour group, SU(2)l x 
SU(2) R , just like the original QCD Lagrangian. 

The trace of the Polyakov loop, $, and its conjugate, will be treated as classical 
field variables throughout this work. In the absence of quarks, we have $ = $ 
and the Polyakov loop serves as an order parameter for deconfinement. The phase 
transition is characterized by the spontaneous breaking of the Z(3) center symmetry 
of QCD. The temperature dependent effective potential U has the following general 
features. At low temperatures, U has a single minimum at $=0, while at high 

^^more precisely: the Polyakov line with periodic boundary conditions 
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temperatures it develops a second one which turns into the absolute minimum above 
a critical temperature T . In the limit T — > oo we have $ — > 1. The function 
£V($,$,T) will be fixed by comparison with pure-gauge Lattice QCD. We choose 
the following general form in accordance with the underlying Z(3) symmetry: 

M (^,p )= _ M r) fa 6, 2 

T 4 2 6 v ; 4 v; v; 



with 



6 2 (T) = a + a 1 (^)+a 2 ^) 2 + a 3 (^) 3 . (6) 



A precision fit of the coefficients a*, 6j is performed to reproduce the lattice data 
(see section 3). 

Using standard bosonization techniques the Lagrangian (1) can be rewritten in 
terms of the auxiliary field variables a and n: 



a +71 



-eff 



-U ($,$,T) -iTi\nS~\ (7) 



2G 

where an irrelevant constant has been dropped and 

S , - 1 = i7^-7o^°-M (8) 

is the inverse quark propagator with 

M = rho — o — 275T • 7? . (9) 

The trace in (7) is taken over colour, flavour and Dirac indices. The field equations 
for a, 7r, $ and $ are then solved in the mean field approximation 2 . The expectation 
value (7?) of the pseudoscalar isotriplet field is equal to zero for isospin-symmetric 
systems. 

The a field has a non-vanishing vacuum expectation value as a consequence 
of spontaneous chiral symmetry breaking. Solving the field equations for a, the 
effective quark mass m is determined by the self-consistent gap equation 

m = m — (a) = m — G^ip). (10) 

Note that (a) = G(^)ip) is negative in our representation, and the chiral (quark) 
condensate is (ipip) = (tp u ipu + 4>dipd)- For later purposes we note that $ and <£> are 
two independent field variables in the general case of finite quark chemical potential 
fi. They become equal in the limiting case fi = 0. Their (thermal) expectation 
values (<£>) and ($) are both real [56] but differ at non-zero //. 

2 In the mean field approximation the fields are replaced by their expectation values for which, 
in later sections, we will continue using the notation a, $ and $ for simplicity and convenience. 
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a 




a 2 


a>3 


h 


h 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 



Table 1: Parameter set used in this work for the Polyakov loop potential (5, 6). 



Before passing to the actual calculations, we summarize basic assumptions be- 
hind eq. (1) and comment on limitations to be kept in mind. In fact the PNJL 
model (1) is quite schematic in several respects. It reduces gluon dynamics to a) 
chiral point couplings between quarks, and b) a simple static background field rep- 
resenting the Polyakov loop. This picture cannot be expected to work beyond a 
limited range of temperatures. At large T, transverse gluons are known to be ther- 
modynamically active degrees of freedom, but they are ignored in the PNJL model. 
To what extent this model can reproduce lattice QCD thermodynamics is nonethe- 
less a relevant question. We can assume that its range of applicability is, roughly, 
T < (2 — 3)T C , based on the conclusion drawn in ref. [57] that transverse gluons 
start to contribute significantly for T > 2.5 T c . 

3 Parameter fixing 

3.1 Polyakov loop effective potential 

The parameters of the Polyakov loop potential U are fitted to reproduce the lat- 
tice data [58] for QCD thermodynamics in the pure gauge sector. Minimizing 
$,T) one has <3> = i> and the pressure of the pure-gauge system is evalu- 
ated as p(T) = —U(T) with $(T) determined at the minimum. The entropy and 
energy density are then obtained by means of the standard thermodynamic rela- 
tions. In Fig. 1 we show the (scaled) pressure, energy density and entropy density 
as functions of temperature. The lattice data are reproduced extremely well using 
the ansatz (5, 6), with parameters summarized in Table 1. The critical temperature 
T for deconfinement appearing in Eq. (6) is fixed at T = 270 MeV in the pure 
gauge sector. The resulting effective potential is displayed in Fig. 2 for two different 
temperatures: T = 200 MeV (below T ) and T = 320 MeV (above T ). 

With the same parametrization, we are also able to reproduce the lattice data [59] 
for the temperature dependence of the Polyakov loop itself. A comparison between 
these data and our results is shown in Fig. 3. The Polyakov loop vanishes below the 
critical temperature T , at which point it jumps discontinuously to a finite value, 
indicating a first order phase transition. It tends to one at large temperatures, as 
expected. 
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Figure 3: Polyakov loop as a function of temperature in the pure gauge sector, compared to 
corresponding lattice results taken from Ref. [59] . 



The pure NJL model part of the Lagrangian (1) has the following parameters: the 
"bare" quark mass m , a three-momentum cutoff A and the coupling strength G. 
We choose to reproduce the known chiral physics in the hadronic sector at T = 
and fix the three parameters by the following conditions: 

• The pion decay constant is reproduced at its empirical value, f n = 92.4 MeV. 
In the NJL model, f n is evaluated using the following equation: 



with the effective (constituent) quark mass m determined self-consistently by 
the gap equation (10). 

• The quark condensate becomes 



3.2 NJL sector 




(11) 



4m 1^ (m) 



(12) 



with 




(13) 



Its "empirical" value derived from QCD sum rules is 



(^> 1/3 - (</W>d} V3 = - (240 ± 20) MeV . 



(14) 
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A [GeV] 


G[GcV- 2 ] 


m [MeV] 


0.651 


10.08 


5.5 


1/3 [MeV] 


/.[MeV] 


m^MeV] 


251 


92.3 


139.3 



Table 2: Parameter set used in this work for the NJL model part of the effective 
Lagrangian (1), and the resulting physical quantities. For these values of the pa- 
rameters we obtain a constituent quark mass m =325 MeV. 



• The current quark mass m is fixed from the Gell-Mann, Oakes, Renner 
(GMOR) relation which is satisfied in the NJL model: 

< - • (15) 

J TV 

In the chiral limit, m = and m n = 0. 

The values of the NJL model parameters, together with the resulting physical 
quantities, are summarized in Table 2. 

4 Results at finite T and fi 

We now extend the model to finite temperature and chemical potentials using the 
Matsubara formalism. We consider the isospin symmetric case, with an equal num- 
ber of u and d quarks (and therefore a single quark chemical potential //). The 
quantity to be minimized at finite temperature is the thermodynamic potential per 
unit volume: 

n (T, //) = U(*,*,T)-*Yf -^Tr In § — M) + ^. 

(16) 

Here u; n = (2n + l)7rT are the Matsubara frequencies for fermions. The inverse 
quark propagator (in Nambu-Gorkov representation) becomes 

cr 1 ( n ° n\ = ( loP° -1 ■p-m--fo(^ + iA 4 ) \ 

{P,P) V lo p°-l-p-m + lo (fi + tA 4 ) J ■ 

(17) 
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Using the identity Tr In (X) = lndet (X) we reduce the trace in (16) and find: 



2 

_ 2N f T j ^{Tr c ln[l + Le-(^/ T ] +Tr c ln[l + Lt e -^ + ^ T ]} 

(18) 

where we have introduced the quark quasiparticle energy E p = \f§ 2 + m 2 . The last 
term involves the NJL three- momentum cutoff A. The second (finite) term does not 
require any cutoff. A small violation of the underlying chiral symmetry at T > 0.4 
GeV, resulting from this procedure, is of no practical relevance since the model is 
supposed to be applied only at temperatures and chemical potential well below A. 
The remaining colour trace is then performed with the result 

In det [1 + L e-( E *-^ T ] + In det [l + iS e~^ + ^ T ] 
= In [1 + 3 ($ + $ e -^-^/ T ) e~^-^' T + e" 3 ^-^] 
+ In [1 + 3 (<5 + $ e -^ + ^/ T ) e - ( ^ + ^ )/T + e - 3( ^ )/T ] . (19) 

From the thermodynamic potential (18) the equations of motion for the mean fields 
a, $ and $ are derived through 

dtt dtt dtt 

*= ' 5f = °- a¥ = ' (20) 

This set of coupled equations is then solved for the fields as functions of temperature 
T and quark chemical potential [i. 

Fig. 4(a) shows the chiral condensate together with the Polyakov loop $ as 
functions of temperature at n — where we find again $ = $. One observes that the 
introduction of quarks coupled to the a and $ fields turns the first-order transition 
seen in pure-gauge Lattice QCD into a continuous crossover. The original 1st order 
transition in the pure-gauge system appears at a critical temperature To = 270 MeV. 
With the introduction of quarks, the crossover transitions for the chiral condensate 
(ipip) and for the Polyakov loop perfectly coincide at a lower critical temperature 
T c ~ 220 MeV (see Fig. 4(b)). We point out that this feature is obtained without 
changing a single parameter with respect to the pure gauge case. The value of the 
critical temperature that we obtain is a little high if compared to the available data 
for two-flavour Lattice QCD [60] which gives T c = (173 ± 8) MeV. On the other 
hand, it is presently being discussed that detailed continuum extrapolation of these 
data can increase this temperature up to 210 MeV [61]. For quantitative comparison 
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Figure 4: Left: scaled chiral condensate and Polyakov loop $(T) as functions of 
temperature at zero chemical potential. Right: plots of d{^)/dT and <9$/<9T. 



with existing lattice results we choose to reduce T c by rescaling the parameter T 
from 270 to 190 MeV. In this case we loose the perfect coincidence of the chiral 
and deconfinement transitions, but they are shifted relative to each other by less 
than 20 MeV. When defining T c in this the average of the two transition 

temperatures we find T c = 180 MeV. This is also consistent with the observations 
reported in [62]. 

As we turn to non-zero chemical potential, we find that $ and $ are different 
from each other, even if they are both real. They will finally coincide again at high 
temperatures, as can be seen in Fig. 5. This feature was already observed in [56]. 

With increasing chemical potential, the crossover pattern evolves to lower tran- 
sition temperatures (see Fig. 6) until it turns to a first order transition around 
fi ~ 0.3 GeV. At this point Cooper pairing of quarks presumably sets in. A more 
detailed discussion of the critical point and its neighbourhood therefore requires the 
additional incorporation of explicit diquark degrees of freedom in the PNJL model. 
Further developments along these lines will be reported elsewhere. 

5 Detailed comparison with Lattice QCD 

A primary aim of this work is to compare predictions of our PNJL model with the 
lattice data available for full QCD thermodynamics at zero and finite //. Consider 
first the pressure of the quark-gluon system at zero chemical potential: 

p (T, ix = 0) = -n (T, fj, = 0; a(T, 0), $(T, 0), 0)) , (21) 
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Figure 5: Averaged sum and difference of $ and $ as functions of the temperature at finite \x. 




Figure 6: Constituent quark mass (a) and Polyakov loop (b) as functions of temperature for 
different values of the chemical potential. 
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where a (T, 0) , $ (T, 0) and $ (T, 0) are the solutions of the field equations at fi- 
nite temperature and zero quark chemical potential. Our results are presented in 
Fig. 7(a) in comparison with corresponding lattice data. We point out that the 
input parameters of the PNJL model have been fixed independently in the pure 
gauge and hadronic sectors, so that our calculated pressure is a prediction of the 
model, without any further tuning of parameters. With this in mind, the agree- 
ment with lattice results is quite satisfactory. One must note that the lattice data 
are grouped in different sets obtained on lattices with temporal extent N t — 4 and 
N t = 6, both of which are not continuum extrapolated. In contrast, our calculation 
should, strictly speaking, be compared to the continuum limit. In order to perform 
meaningful comparisons, the pressure is divided by its asymptotic high-temperature 
(Stefan-Boltzmann) limit for each given case. At high temperatures our predicted 
curve should be located closer to the N t = 6 set than to the one with Nt = 4. 
This is indeed the case. Furthermore, Fig. 7(b) shows the predicted "interaction 
measure", (e — 3p)/T 4 , in comparison with lattice data for N t = 6. One should of 
course note that the lattice results have been produced using relatively large quark 
masses, with pseudoscalar-to- vector mass ratios mps/my around 0.7, whereas our 
calculation is performed with light quark masses corresponding to the physical pion 
mass. We have investigated the dependence of the pressure and of the energy density 
on the quark mass and found that the critical temperature scales approximately as 
T c ~ Tjyjn^ = 0) + 0.04 m^, in agreement with the behaviour found in [60]. Once 
the rescaling of T c is taken into account, the curves plotted in Figs. 7(a) and 7(b) 
as functions of T/T c have negligible remaining dependence on the quark mass. 

At non-zero chemical potential, quantities of interest that have become accessible 
in Lattice QCD are the "pressure difference" and the quark number density. The 
(scaled) pressure difference is defined as: 

AP(3» = p(T,ti)-p(T,ti = 0) 

A comparison of Ap, calculated in the PNJL model, with lattice results is presented 
in Fig. 8. This figure shows the scaled pressure difference as a function of the 
temperature for a series of chemical potentials, with values ranging between /i = 
0.2 Tc°^ and /i ~ Tc°\ The agreement between our results and the lattice data is 
quite satisfactory. 

A related quantity for which lattice results at finite \i exist, is the scaled quark 
number density, defined as: 

n q (r,/z) = i an(r,/z) 

T 3 T 3 d/ji 1 ' 

Our results for function of the temperature, for different values of the quark 

chemical potential, are shown in Fig. 9 in comparison with corresponding lattice data 
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Figure 7: (a) Scaled pressure divided by the Stefan-Boltzmann (ideal gas) limit as a function 
of temperature at zero chemical potential: comparison between our PNJL model prediction and 
lattice results corresponding to N t = 4 and N t = 6. (b) Scaled interaction measure compared to 
lattice results for N t — 6. In both cases, the lattice data are taken from Ref. [63] 
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Figure 9: Scaled quark number densities as a function of temperature at different values of the 
chemical potential, compared to lattice data taken from Ref. [4]. 

[4]. Also in this case, the agreement between our PNJL model and the corresponding 
lattice data is surprisingly good. 



It is instructive to study the effect of the Polyakov loop dynamics on the behaviour 
of the quark density n q . The coupling of the quark quasiparticles to the field <3> re- 
duces their weight as thermodynamically active degrees of freedom when the critical 
temperature T c is approached from above. At T c the value of $ tends to zero and 
the quasiparticle exponentials exp[— (E p ± fi)/T] are progressively suppressed in the 
thermodynamic potential as T — > T c . This is what can be interpreted as the impact 
of confinement in the context of the PNJL model. In contrast, the standard NJL 
model without coupling to the Polyakov loop does not have this important feature, 
so that the quark density leaks strongly into the "forbidden" domain T < T c ~ 170 
MeV, as demonstrated in Fig. 10. 

It is a remarkable feature that the quark densities and the pressure difference at 
finite /i are so well reproduced even though the lattice "data" have been obtained by 
a Taylor expansion up to fourth order in fi, whereas our thermodynamic potential is 
used with its full functional dependence on fi. We have examined the convergence in 
powers of fi by expanding Eq.(18). It turns out that the Taylor expansion to order 
fi 2 deviates from the full result by less than 10 % even at a chemical potential as 
large as fi ~ T c . When expanded to 0(fi A ), no visible difference is left between the 
approximate and full calculations for all cases shown in Figs. 8 and 9. 
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T/T 



Figure 10: Comparison between the results in the PNJL model (solid line) and in the standard 
NJL model (dashed line) for the quark number density at \i = 0.6T C . The effect of the missing 
confinement is evident in the standard NJL model. 

6 Summary and conclusions 

We have studied a Polyakov-loop-extended Nambu and Jona-Lasinio (PNJL) model 
with the aim of exploring whether such an approach can catch essential features of 
QCD thermodynamics when confronted with results of lattice computations at finite 
temperature and non-zero quark chemical potential. This PNJL model represents a 
minimal synthesis of the two basic principles that govern QCD at low temperatures: 
spontaneous chiral symmetry breaking and confinement. The respective order pa- 
rameters (the chiral quark condensate and the Polyakov loop) are given the meaning 
of collective degrees of freedom. Quarks couple to these collective fields according 
to the symmetry rules dictated by QCD itself. 

Once a limited set of input parameters is fitted to Lattice QCD in the pure 
gauge sector and to pion properties in the hadron sector, the quark-gluon thermo- 
dynamics above T c up to about twice the critical temperature is well reproduced, 
including quark densities up to chemical potentials of about 0.2 GeV. In particular, 
the PNJL model correctly describes the step from the first-order deconfinement tran- 
sition observed in pure-gauge Lattice QCD (with T c ~ 270 MeV) to the crossover 
transition (with T c around 200 MeV) when Nf — 2 light quark flavours are added. 
The non-trivial result is that the crossovers for chiral symmetry restoration and 
deconfinement almost coincide, as found in lattice simulations. The model also re- 
produces the quark number densities at various chemical potentials remarkably well 
when confronted with corresponding lattice data. Considering that the lattice re- 
sults have been found by a Taylor expansion in powers of the chemical potential, 
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this excellent agreement came as a surprise and indicates rapid convergence of the 
power series in //. 

Further developments will be directed towards improvements to overcome some 
obvious limitations. First, the NJL model operates with a constant four-point cou- 
pling strength which supposedly averages the relevant running coupling over a lim- 
ited low-energy kinematic domain, corresponding to temperatures T < 2T C and 
chemical potentials fi < 0.3 GeV. Contacts with the high-temperature limit of QCD 
and the HTL approaches need to be established. Secondly, in order to proceed into 
the range of larger chemical potentials, diquark degrees of freedom need to be ex- 
plicitly involved. Also, the effective potential for the Polyakov loop field, determined 
so far entirely as a function of temperature by investigating the pure gauge sector, 
must be examined with respect to its dependence on the chemical potential. And 
furthermore, the extension to 2+1 flavours with inclusion of strange quarks must be 
explored. 

Nevertheless, considering the simplicity of the PNJL model, the conclusion that 
can be drawn at this point is promising: it appears that a relatively straightforward 
quasiparticle approach, with its dynamics rooted in spontaneous chiral symmetry 
breaking and confinement and with parameters controlled by a few known properties 
of the gluonic and hadronic sectors of the QCD phase diagram, can account for 
essential observations from two-flavour N c = 3 Lattice QCD thermodynamics. 
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